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Abstract 

It is shown that in equilibrium a canonical ensemble of particles with two- 
particle interaction the Gibbs distribution function may be expressed uniquely 
through a pair distribution function. It means, that for given values of the par- 
ticle number N, volume V, and temperature T, the pair distribution function 
contains as many information about the system as a full Gibbs distribution. 
The latter is represented as a series expansion in the pair distribution func- 
tion. A recurrence relation system is constructed, which allows all terms of 
this expansion to be calculated successively. 



1 Introduction 

In classical statistical mechanics all properties of a closed system of N particles in 
a volume V (canonical ensemble) are described by a distribution function Dn{X), 
where X is a set of system phase variables. For a system with an additive particle 
interaction, multi-particle distribution functions are introduced [TJ Ej to calculate 
thermodynamical values of the system. These multi-particle functions are consid- 
ered to have less information than the full Gibbs iV-particle distribution function. 
Thereby the smaller the order of the distribution function is the less information it 
contains 0. But there does not exist a proof of this statement in the literature. 

At the same time it is known that for an equilibrium system of A" non-interacting 
particles in the volume V with the temperature T the GDF decomposes into a 
product of one-particle distribution functions. This means that all information about 
the system is contained in the one-particle distribution function. 

In this article a problem on equivalence of describing an equilibrium system by 
full Gibbs distribution function (GDF) and pair (two-particle) distribution function 
(PDF) is considered. It is used a non-normalized GDF and particle distribution 
functions are defined trough it. The motive for that will be explained in conclud- 
ing section. As to a description of system the non-normalized GDF contains all 
information about it. 
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We shall show that for the equilibrium system of iV pair-interacting particles in 
the absence of an external field, when the energy Un equals 

£Mqi,---,qjv) = M 2(qi,qj), (l) 

l<i<j<N 

a Gibbs distribution may be expressed uniquely through a PDF. It means the PDF 
contains all information on the system under consideration. 



2 Formulation of the problem 

The equilibrium canonical ensemble is described by the Gibbs distribution function 1 

D N (q 1 ,...,q N ) = exp{- / g[/ A r(qi,---,qjv)}, G V, i = l,...,N, ft = 1/kT. 

(2) 

where k is the Boltzmann constant, V is a domain of three-dimensional Euclidean 
space 2 . Here and below in this paper we shall suppose that the configuration vari- 
ables q belong to the domain V. 

For the system of particles with the energy (PJ the expression (J2J) may be rewrit- 
ten as follows 

D N (q_i,...,ci N ) = Yl u (<li>qi)> ( 3 ) 

l<i<j<N 

where we introduce the function 

V (q,q') d =exp{-/3 M2 (q,')}. (4) 

Due to reality of the potential «2 this function is nonnegative. And by virtue of 
symmetry with respect to its arguments the function f(q, q') is symmetrical too. 
It follows from the expression (j3J) that statistical properties of the system under 
consideration are described completely by specifying a single function v (q, q') of 
two configuration variables. Naturally, there appears an assumption that, if this 
system may be described entirely with the help of one function of two variables 
f(q, q'), then it may also be described with the help of another function of two 
variables connected with previous one, namely a pair distribution function. We 
prove this assumption here. 

Define s-particle distribution functions Q s starting from the non-normalized GDF 

© 

&(qi, . . . , q s ) = ( N _ s yy N J D n(<Ii, ■ ■ ■ > qv)dq s +i • • • dq N s = l,...N. (5) 



^^Here we consider only a configuration part of the distribution function, because its momentum 
part is factorized per a one-particle Maxwell function. 

2 Three-dimensionality of space is not essential. It may be of any finite dimensionality. Only 
fmiteness of the Lebegue measure (volume) of the manifold V is important. 
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Here V is a volume of the domain V. The integration will be carried out over the V 
and the domain of integration will not be indicated explicitly. 

Due to the properties of v (q, q') all multi-particle distribution functions are real, 
nonnegative, and symmetrical. These functions are related to the ones f s , introduced 
by Balescu [I], and to F s introduced by Bogoliubov [2] as follows: f s = Q~^V N Q S and 
F s = (Q N N\)-\N - s)\V N - s g s respectively, where Q N is a configuration integral 

f V N (N-s)\ f 
Q N = dq 1 --- dq N D N ( 1: ...,q N ) = — / dq x • • • dq s Q s ( u . . . , q s ). (6) 

Write down explicitly an expression for the PDF 

&(qi,q2) = N(y ^N - I II v(qi,qj)dq 3 - ■ -dq N . (7) 

l<i<j<N 

Let us introduce a nonlinear operator $ as follows 

= yt^2 / n v ^ ^^ q3 ■ ■ ■ dqN - ( 8 ) 

l<i<j<N 

Then the relation (jJJ) may be considered as a nonlinear operator equation relative 
to v 

-<S>(v) d ^F(v,g 2 ) = 0. (9) 
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N(N-l) 
If this equation has a solution 

v = X(S 2 ), (10) 

then its substitution into (j2J) permits one to express the non-normalized Gibbs dis- 
tribution Z)jv(qi, • • • , qw) through the two-particle function ^(q, q')- Accordingly 
all s-particle functions may also be expressed via ^(q, q')- 

For calculating the function v from equation © it is necessary to set an addi- 
tional condition 

where and v ^ are known functions. It is easy to find such a condition assuming 
particularly 

t,<°W) = i, ^ 0) (q,q') = ^^, (12) 

which corresponds to the choice of noninteracting particle system as a reference 
one 3 . 



3 Evidently, instead of unity for (q, q') we may take an arbitrary constant C. In this case 
an expression for {<t\, q') is calculated easily too. 
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3 Proof of the existence of a solution 



In functional analysis there exist a number of theorems about implicit functions for 
operators of different smoothness degree (e.g. [HIE]). We use a theorem for analytical 
operators in Banach spaces in the form presented in the book [2 Theorem 22.2]. 
Theorem. Let T(v,Q 2 ) be an analytical operator in D r (v^°\ Ei)+D p (Q^\ E) 

with values in E 2 - Let the operator B = f —dJ-{v^°\ Q^)/dv have a bounded inverse 
operator. Then there exist such positive numbers p\ and r\ that in the solid sphere 
D n (v(°\Ei) there exists a unique solution \1(J\) to equation (0J). This solution is 
defined in the solid sphere D Pl (Q^\ E) , being analytical in it and satisfying the 
condition = xiS^)- 

Here D r (xo,£) is a solid sphere of the radius r in the vicinity of an element xo 
of Banach space S, the sign "+" denotes an algebraic sum of manifolds, and dT /dv 
is a Frechet derivative |H Ej of a nonlinear operator T . 

To investigate the problem of existence and uniqueness of the solution of equation 
Q it is necessary to verify fulfillment of the conditions of the theorem for the system 
under consideration. 

At first define spaces E\ and E, in which the operator T is set, and space E 2 
containing a range of values of this operator. We shall use a class of functions v (q, q') 
bounded almost everywhere in V. The space of such functions becomes a complete 
normalized (Banach) space of essentially bounded functions L OT if we introduce the 
following norm [H E] 

||f|| = vraisup \v(q, q')|, (13) 

(q,q')eVxV 

where "vraisup" denotes an essential exact supremum of the function on the indi- 
cated set, and VxV is a direct product of the manifold V by itself, that is a set of 
ordered pairs (q, q'). Moreover f(q, q') is a symmetrical function of its arguments. 
Therefore v (q, q') G (VxV). The index (s) means that we use a space of symmet- 
rical functions. Such properties of the function t>(q, q') defined by (J3J) are provided 
by any interaction potential bounded below almost everywhere in the domain V. 
All potentials conventionally used in statistical mechanics satisfy this condition. 

Using the properties of real functions and their integrals one can show that 
for any functions v from (VxV) both GDF and all multi-particle ones are also 
essentially bounded functions being defined on appropriate sets. Thus the space of 
functions v is E\ — L^(VxV), the space of functions Q 2 is E = L^(VxV), and the 
space containing the range of values of the operator T is E 2 = L^(VxV). The set 
of pairs (v , Q 2 ) of nonnegative functions in the direct product E\ x E is a range of 
definition 4 of the operator J 7 , with its range of values belonging to E 2 . Everywhere 
in the range of definition the operator T{v , Q 2 ) is analytical. 

4 In general, the requirement on nonnegativity of the functions follows from physical reasons. 
From the mathematical point of view the operator J-(v 1 Q2) © is defined on a whole space E\ x E. 
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Calculate an operator B defined in the theorem formulation for chosen values of 
the functions v {0) and (321). It is a linear operator acting from space E\ to space 
E 2 . Using the first equality in Q, we obtain the expression 



B 



d®(v 



{0)\ 



dv 



(14) 



Substitute into it the expression (JHJ) for the operator $(u) and use the condition 
v(°> = 1. We obtain the operator B as follows 

(Bh)(q, q') = h(q, q') + (N - 2) p(q) + fc(q')] + ~ 2 X^ " 3 ) ft, (15) 

where we introduce designations 

Mq) = ^J h(q,q')dq>, h = ^j h{q^)dqdq'. (16) 

In these expressions h e 5/i G i?2- The range of definition of B is a whole space 
£?x : D(B) = Ei and its range of values is a space E 2 : R(B) = E 2 . The operator B 
is bounded. Its norm defined by (see eg. [HIS]) 



\B\\ = sup \\Bh\ 

{heEi;\\h\\=l} 



(17) 



is evaluated easily with the help of (|15|). Using properties of norms we obtain the 
following inequality 



\Bh\\<\\h(q,q')\\ + (N-2)\\h(q) + h(q')\\ + 



(N - 2)(N - 3) ..= 



;is) 



Evaluate the norms in right hand side of (fTB|) . The value /i is a constant. For its 
norm we have an estimate 



\h\ 



1 



h(q, q')dqdq' 



1 

V 1 



< 



1 



vrai sup \h(q, q') \ dqdq' 

(q,q')GVxV 



I /i || dqdq' = \\h\ 



(19) 



For the norm of a sum of the functions h(q) and h(q') we obtain the following 
estimate 

\\h(q) + h(q')\\ = vrai sup \h(q) + h(q')\ < 2 vrai sup p(q) = 

(q,q')£VxV qGV 



— vrai sup 

V qGV 



h{q, p)dp 



< 



V 



vrai sup \h(q, p)| dp < 

qGV 
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2 r 

< — / vrai sup \h(q, p)| dp — 2\\h\ 

V J ( q ,p)gVxV 



(20) 



Substituting (JT§j) and (J2Uj) in the inequality (JTHJ), we obtain the following estimate 
for the norm of the function (.B/i)(q, q') 

\\Bh\\< N{N ~ % h\\. (21) 

Taking into account the definition (|T7jl . we obtain an estimate of the norm of the 
operator £> 

„ „ N(N-l) , s 

LB < — 1 -. 22 

II II - 2 \ ) 

One can show that actually in (|22j) the sign of equality is valid. Indeed because 
there is at least one function h(q, q') G E\ with the norm equal to unity for which 
\\Bh\\ equals the value in the right hand side of (}2"Tj) . namely h(q, q') = 1, then the 
inequality (|2*2*j) reduces to the equality 

1*1 = SiSfSL. (23) 

Now define a kernel (a space of zeros) of the operator B, i.e. find a solution to 
the homogeneous equation Bh = 0. Using (JT3J), we obtain the equation 

Kq, q) + (N- 2)[h(q) + fc(q')] + ~(iV - 2)(iV - 3)S = 0. (24) 

Integration of this equation over q' and then over q allows one to obtain its solution. 
It turns out to be trivial h(q, q') = 0. 

This means that the corresponding inhomogeneous equation Bh = g is uniquely 
solvable. Its solution has the form 

Mq,q') d = f (B-^)(q,q') = <?(q, q') - ^|[^(q) + ^(q')] + ^j^l (25) 

where g(q) and ^ are determined by the formulae analogous to (JTHj) where /i is 
replaced by g. It is easy to see that the operator B~ x is specified for all g G E 2 , that 
is indeed R(B) = E 2 . It is easy to verify that the conditions B~ l B = BB^ 1 = I 
are valid. In other words, the operator B is both left and right inverse. Therefore 
it is the inverse operator to B j^j. 

The operator B~ x is bounded. For its norm an upper bound is easily evaluated 
just as for the operator B. At once write down this estimate 

i 3iV 2 — 6iV + 2 
LB < ; i—. 26 

Thus all conditions of the theorem are fulfilled. Hence there is a unique solution 
(I10J) expressing unambiguously the function u(q, q') in terms of PDF Qi{*\.-, q')- This 
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solution is determined in the solid sphere D Pl (Qf\ E), analytical in it, and satisfies 
the conditions (fTTj) . (fT2*|) . We shall not produce here a proof of this theorem. It is 
presented in many books in functional analysis. Indicate merely the procedure of 
constructing a solution to equation (jUJ) and estimating the convergence domain size. 
Introduce functions g(q, q') and h(q, q') by the relations 

g 2 = N{N v ~ 1 \ l + g), v = l + h. (27) 

Here and below /i(q, q') is a Mayer function. Equation Q with the operator $(t>) 
(JB1) may be rewritten as 

N 

h = B- l g + B- lS ^B s (h), (28) 

s=2 

where M = N(N — l)/2 and B s (h) is a homogeneous power operator of the order s 
with respect to h equal to 

B sW = ~-^ [ dq 3 ---dq N h{X Kl )...h(X Ks ). (29) 

•* l<Kx<-<Ks<ST 

Here we introduce a generalized index K varying from 1 to M enumerating all 
possible different index pairs (ij), satisfying the condition % < j which can be 
constructed from the numbers 1, 2, . . . , N. Such one-to-one correspondence K 
(ij) is easily found. Xk designates the corresponding pair of coordinates (q^q,). 
One can consider that the power operator B s (h) is generated by the s-linear operator 
R s (h u ...,h s ) El 

R s (h, ...,h s ) = - y\ N j dq 3 --- dq N ^ h x (X Kl ) ■ ■ -h 3 (X Ka ). (30) 

l<Ki<-<K s <N 

Using this operator, we rewrite equation ()28p as follows 

Af 

h = B~ l g + B- lS ^R s (h } ...,h). (31) 
Its solution is found as a power series in g 

oo 

h = J2h(g), (32) 



k=l 



where hk(g) is an homogeneous power operator of the order k acting from space 
E 2 to Ei. This means that each operator hk(g) is a function of two configuration 
variables. 
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Substituting the expression (J32|) into equation (f3~T|) and producing some rear- 
rangement of terms in the right-hand side, we derive the following relationship 



J2hk(g) = B- 1 g + B- 1 ^^ 9 ^ E Rs(hkM,---,h hs (g)), (33) 

fc=l k=2 s=2 fciH \-k s =k 



where 



1, s<Af, 
0, s > Af. 



In accordance with the theorem on uniqueness of analytical operators [7j, equating 
terms of identical orders on g, we obtain a recurrent system for determining the 
operators hk{g) 

h 1 (g) = B- 1 g, 

k 

h k (g) = B- 1 J2 sX Yl Rs(h kl ,...,h ks ), k = 2,3,.... (34) 

s=2 fciH \-k s =k 

All expansion terms of (|3*2|) are calculated successively from this system. And thus 
a formal solution h(g) is derived in the form of the above series (}3*2|) . It is also 
necessary to prove a convergence of this series. It is established, just as in [3], by 
Goursat method [S] of majorizing functions. 
Consider a subsidiary equation 



f N 

l s=2 



A! (35) 



s!(JV-s)! 



in which we designate £ = \\h\\, t] = \\g\\ and 7 = ||-B _1 ||. The right-hand side of this 
equation majorizes the norm of the right-hand side of (|3ip. A solution to equation 
()35j) can be found in the form of expansion in series with respect to 77: £ = Y1T=1 £kV k - 
For coefficients £fc it is easy to derive a recurrent system analogous to the structure of 
(jnij). Just as in a classical case [S] it is easy to show that ||^fc(fi , )|| < £,kV h - Therefore 
if the series for £ converges, then the series (}3*2|) absolutely converges. 
On the other hand, equation (J3*5j) can be rewritten as 

7 T7 = £- 7 [(l + A/ '-i-A/'e (36) 

We cannot write out a solution to this equation explicitly by quadratures, but it 
is easy to show graphically that the solution £(77) satisfying the condition £(0) = 
exists in the region rj < rj . This solution is unique, analytical and monotonically 
increasing. The value £0 determines the maximum point of the function in the right- 
hand side of equation (|3l)|) and equals £0 = (1 + l/jN) 71 - 1 ■ The value rj is defined 
by this function maximum and equals 



% = (JV-1)(1+ 1 



7AA 



1 



7-AT 
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Thus the series (jH2j) converges if \\g\\ < r/ . Thereby, as a radius p\ of the solid 
sphere from the theorem, any number smaller than 7/0 may be chosen. By virtue 
of a monotonic increase of the function £(77), the radius 7*1 may be taken equal to 
7*1 = £(pi). So the solution (fTUj) to the problem Q exists, it is unique and may 
be presented in the form of the expansion (|3"2)l into an absolutely convergent series 
with respect to g — (Q 2 — G2 )/&2 a ^ l eas t in the region \\g\\ < r/ . 

Since the Gibbs distribution is a polynomial of the order M relative to the 
function h, then it is analytical function of g and may be expanded in a power series 
with respect to g 

00 

D N = Y,d k (g), (37) 

where dk(g) is a homogeneous power operator of the order k acting from a space 
E 2 into a space En = £00 (Vx • • • xV), that is, each operator dk(g) is a function 



N co-factors 



of N configuration variables. Due to determining the PDF (J7|) and the theorem on 
uniqueness of analytical operators [7] these functions satisfy the conditions 

do = 1, 

1 



yN- 



dq 3 ■ ■ ■ dqjvdi(qi, . . . , q^) = g(qi), (3J 



j dqt ■ ■ ■ dq N d k (qi, . . . , q^) = 0, k = 2,3, 

Substituting the expansion (jS2j) into into (jSJ) and producing some rearrangement 
of terms, we can derive uniquely all terms d k {g) of the expansion ([HZ} in the form 

k 

d k (g) = Y, 9 ^ Yl h kl (X Kl )---h ks (X Ks ), k — 1,2,..., 

s=l fcjH h<fe a =fc l<E'i<---<A' a <A^ 

(39) 

where hi(XK.) are expressed in terms of g by formulae ([34)1 . The expressions ob- 
tained determine an expansion of the GDF in series with respect to PDF. 



4 Concluding remarks 

Using the theorem of existence and uniqueness of the solution to equation ©, we 
reveal that Mayer function h(q, q') may be expressed uniquely through the PDF 
^2(q ; q') ©• It permits us to present the non-normalized Gibbs distribution as a 
nonlinear operator function of PDF and write its expansion to power series. The 
procedure of successive calculation of the all terms of this series is formulated. 

So knowing a pair distribution function, we can express the Gibbs distribution 
uniquely in terms of this function. In turn PDF is expressed uniquely via GDF 
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(see (jSJ). This means that descriptions of equilibrium canonical ensemble by both 
GDF and PDF are equivalent, at least for h belonging to the solid sphere D ri (0, E\). 
Hence PDF contains the same quantity of information as the full GDF, as well as all 
s-particle distribution functions for s > 2 which are also expressed uniquely through 
PDF. 

The one-particle distribution function £/i(q) does not seem to contain all infor- 
mation about the system with two-particle interaction. If one tries to realize the 
above programme, then one fails to satisfy all conditions of the theorem on existence 
and uniqueness of the solution. The kernel (a space of zeros) of the corresponding 
operator B consists of arbitrary functions h(q, q') G , (V x V) limited only by one 
condition h(q) = 0. Therefore it is possible for such operator B to have no inverse 
operator. And indeed it has not. 

Note here that for the normalized GDF an analogous programme cannot be 
realized too. In this case the corresponding operator B has also a nontrivial kernel, 
namely a set of all functions constant on V. Mathematically, it is clear from the form 
of the normalized GDF in which the function v may be multiplied by an arbitrary 
constant without changing the distribution function. Physically, it corresponds to 
the fact that the interacting potential is defined with an accuracy to an arbitrary 
constant. 
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